home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / dsp / dr.bub / 96000.lha / 96000 / appb / b110.asm < prev    next >
Assembly Source File  |  1992-04-28  |  3KB  |  49 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.10    N Complex Updates  
  8. ;   dr(I)+jdi(I) = {cr(I)+jci(I)}+{ar(I)+jai(I)}*{br(I)+jbi(I)}, I=1,...,N 
  9. ; 
  10. ;   dr(I) = cr(I) + ar(I) * br(I) - ai(I) * bi(I) 
  11. ;   di(I) = ci(I) + ar(I) * bi(I) + ai(I) * br(I) 
  12. ; 
  13. ;       D5 = ar    D4 = ai    D6 = br    D7 = bi 
  14. ; 
  15. ;       X Memory Organization            Y Memory Organization 
  16. ;                  .                                        .      
  17. ;                 ci2                                      di2     
  18. ;                 cr2                                      dr2     
  19. ;                 ci1                                      di1     
  20. ;        R1 ?     cr1     CADDR    R5 ?     dr1         DADDR 
  21. ;                  .                                        .      
  22. ;                  .                                        .      
  23. ;                 ai2                                      bi2     
  24. ;                 ar2                                      br2     
  25. ;        R0 ?    ai1                                      bi1     
  26. ;                 ar1         AADDR    R4 ?    br1     BADDR 
  27. ;                                                                         Program    ICycles
  28. ;                                                           Words 
  29.      move   #aaddr+1,r0                                       ;    1          1 
  30.      move   #3,n0                                             ;    1          1 
  31.      move   #baddr,r4                                         ;    1          1 
  32.      move   #caddr,r1                                         ;    1          1 
  33.      move   #daddr-1,r5                                       ;    1          1 
  34.      move                        x:(r0)-,d4.s   y:(r4)+,d6.s  ;    1          1 
  35.      fclr   d2                   x:(r0)+n0,d5.s y:(r5),d0.s   ;    1          1 
  36.      do #n,end                                                ; 2          3 
  37.      fmpy d5,d6,d2 fadd.s d2,d0  x:(r1)+,d1.s   y:(r4)+,d7.s  ;    1          1 
  38.      fmpy d4,d7,d2 fadd.s d2,d1  x:(r1)+,d0.s   d0.s,y:(r5)+  ;    1          1 
  39.      fmpy d4,d6,d2 fsub.s d2,d1  x:(r0)-,d4.s   y:(r4)+,d6.s  ;    1          1 
  40.      fmpy d5,d7,d2 fadd.s d2,d0  x:(r0)+n0,d5.s d1.s,y:(r5)+  ;    1          1 
  41. end
  42.                    fadd.s d2,d0                               ;    1         1 
  43.      move                                       d0.s,y:(r5)+  ;    1          1 
  44. ;                                                            -------------------
  45. ; 
  46. ;                                                        Totals: 15        4N+12
  47. ;                                                                (13        4N+10)
  48.  
  49.